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We present a combined study of the angle-resolved-photoemission 
spectroscopy (ARPES) and quantum Monte Carlo simulations to 
propose a novel polaronic metallic state in underdoped cuprates. 
An approximation scheme is proposed to represent underdoped 
cuprates away from 1/2 filling, replacing the many-body Hamil- 
tonian by that of a single polaron with effective electron-phonon 
interaction (EPI), that successfully explains many puzzles such 
as a large momentum-dependent dichotomy between nodal and 
anti-nodal directions, and an unconventional doping dependence 
of ARPES in the underdoped region. 



PACS numbers: 71.10.Fd,71.38.-k,02.70.Ss 



It is established that the physics of high temperature 
superconductors is that of hole doping a Mott insulator 
[H, 0] where even a single hole is substantially influenced 
by many-body effects [3j . Understanding of the dynamics 
of holes in Mott insulators has attracted a great deal of 
interest [H, [H H|. The major interactions are electron- 
electron interactions (EEI) and EPI. The importance of 
the former is no doubt essential since the Mott insulator 
is driven by this EEI, while the latter was considered 
to be largely irrelevant to superconductivity based on 
the observations of a small isotope effect on the optimal 
T c [B[ and an absence of a phonon contribution to the 
resistivity @. 

On the other hand, there is now accumulating evidence 
that the EPI plays an important role in the physics of 



cuprates [lHH- In particular, EPI manifests itself in (i) an 
isotope effect on superfluid density p s and T c away from 
optimal doping Q, and (ii) neutron and Raman scat- 
ments showing strong phonon soft- 



termg 

ening with both temperature and hole doping, indicating 
that EPI is strong [13[. Furthermore, the recent advance- 
ment in the energy /momentum resolution [3| of ARPES 
resulted in the discovery of the dispersion "kinks" at 
around 40-70meV measured from the Fermi energy, in 
the range of the relevant oxygen related phonons 15- 
These particular phonons - oxygen buckling and 
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half-breathing modes are known to soften with doping 
HI El and with temperature [IMl EMI and there 



is mounting evidence relating the kink to the phonon 



anomaly of the bond stretching phonons 
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quick change of the velocity can be predicted by any in- 
teraction of a quasiparticle with a bosonic mode, either 
with a phonon because of EPI 16, 17, 20, |21| or with a 
collective magnetic resonance mode 22J, |23[ . Early stud- 
ies of the doping dependence of the kink revealed a "uni- 
versality" of the kink energy for La2- ;c Sr a ;Cu04 (LSCO) 
over the entire doping range |24| and casted doubts on 
the validity of the latter scenario as the energy scale of 
the magnetic excitation changes strongly with doping. 
More recent studies showed that there is a subtle doping 
dependence of the kink energy |25| which, nevertheless, 
can be explained within the framework of EPI scenario 
as well [21 1 . 

As a more direct test distinguishing between these 
two scenarios, an observation of an isotope effect just 
in the vicinity of the kink 26j and its explanation in 
terms of EPI 27J has given an additional argument in 
favor of a phononic origin of the kink. Furthermore, it 
was shown that the kink is observed in electron-doped 
cuprates around 40-70meV, i.e. at the same energies as 
the hole-doped ones 28] . Hence, it is clear that the kink 
is due to EPI with phonon modes, which are roughly at 
the same energies in hole- and electron-doped cuprates, 
and not because of coupling to a magnetic mode whose 
energy in the electron-doped compounds is not larger 
than 10 meV [H,^. 

Generally, there are two possible pictures of the EPI. 
One is the Migdal-Eliashberg picture where the kinetic 
(Fermi) energy of the electrons ef is much larger than 
the phonon energy 57. In this limit, the multi-phonon 
processes represented by the vertex correction are sup- 
pressed by the Fermi degeneracy and reduced by the adi- 
abatic factor fl/cp, i.e., EPI is basically in the weak cou- 
pling region j3ll [32| ] . The other limit is the polaron picture 
where Ef is much smaller than VL where the multi-phonon 
processes can, in principle, lead to the small polaron for- 
mation. The latter picture of the strong coupling limit 
of a polaron in undoped materials, where ARPES corre- 
sponds to the single hole dynamics in a Mott insulator, 
has been established by a detailed comparison between 
experiment f33T - l35j and theory 36-3!|. The picture ob- 
tained there is that the "quasi-particle" peak observed 
experimentally is that of Franck-Condon multiphonon 
band, which follows the energy dispersion of bare t-J 
model without EPI, while the zero-phonon line has only 
a very small weight Z << 1. Namely the hole polaron is 
in the strong coupling small polaronic state. 

At very small dopings it is clear that the Fermi energy 
Ef of a few holes doped into the Mott insulator is smaller 
than the relevant phonon energy 51. Then the polaronic 
energy dispersion can be inverted and the ARPES disper- 
sion is cut-off by the (small) Fermi energy e p (measured 
from the hole side). In this case the adiabatic factor ft/sp 
is not small and one cannot rely on the weak coupling 
Migdal-Eliashberg approach and, hence, has to treat the 
EPI in the polaron limit including all vertex corrections 



and without considerable approximations. The polaron 
paradigm is complementary to the conventional Migdal- 
Eliashberg picture for metals 31 , 32| , the latter of which 
fails to explain various anomalous features in the under- 
doped region. Therefore it is reasonable to apply the 
polaron picture in the underdoped region. 

In the present paper we interpret the kink in under- 
doped cuprates as a result of a short-range Holstein-like 
interaction of a hole with optical phonons. Since the po- 
laron scenario is more appropriate in the underdoped re- 
gion than the Migdal-Eliashberg one, we approximate the 
interacting many-body Hamiltonian by the single hole 
polaronic Hamiltonian with an effective dimensionless 
coupling constant A of EPI. The effective coupling con- 
stant A in the weakly doped compound is renormalized 
in comparison with the bare constant Ao, corresponding 
to the case of a single hole, by EEL It is known for gas of 
Frohlich polarons both theoretically and experimentally 
ioT - 44 1 that the basic properties of polarons behave with 
increasing concentration as if the coupling constant of the 
EPI is effectively reduced. Comparison with experiment 
shows that effects of weak EPI and Coulomb repulsion 
can be factorized in 2D systems in the light doping limit 
and, besides, this factorization is approximately valid up 



to a doping level i » 0.1 [40j, |4l|. Since the radius of 
Holstein polarons is smaller than that of Frohlich ones 
we expect in our case better applicability of the model- 
ing of the influence of doping by changing the effective 
EPI constant A. 

The above assumption indicates that A should be re- 
garded as the effective coupling constant in the trial 
Hamiltonian which best mimics the experimental results. 
Such an approach has been already applied to the expla- 
nation of the doping dependence of the optical absorp- 
tion spectra of underdoped cuprates and it was shown 
that the effective coupling constant A decreases with in- 
creasing doping [Hj]. 

Using the approximation-free Diagrammatic Monte 
Carlo (DMC) method [H, Sfil we have calculated 
the momentum dependence of the Lehman spectral func- 
tions of a single hole in the extended t — J model, i.e., 
t-t' -t" - J -ph model (Fig. [TJ. Within this model a 
single hole in an antiferromagnet can hop up to the 3rd 
nearest neighbor (t(l) = t, t(2) = t', and t(3) = t")\ 
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Here Cj a is a projected (to avoid double occupancy) 
fermion annihilation operator, rij = ^ a c\ a Ci a < 2 is 
the occupation number, is spin-1/2 operator, (ij,n) 
denotes neighboring sites on n-th coordination sphere 
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(circle) in a two-dimensional lattice. We note, that al- 
though long-range antiferromagnetic (AFM) order and 
AFM gap is destroyed at rather low doping levels, over- 
damped spin waves and short range AFM correlations, 
showing the largest interpla y w ith short-range EPI, per- 
sist up to fairly high doping |50l l5l| . Hence, the most im- 
portant part of the interpaly between short-range AFM 
and short-range EPI survives at moderate dopings. 

The hole also interacts with dispersionless optical 
phonons ft via a short range Holstein coupling 7 52, 53| 



EPI 



H 



- 1/2 7E 

k,q 



k 



[fihu. 



(2) 



(Here 6k and /ik are phonon and hole annihilation op- 
erators, respectively.) We define the dimensionless cou- 
pling constant A = ^ /Atft of the EPI. Our calculations 
are done for parameters which are required to reproduce 
the observed dispersions 0H3 in Sr 2 Cu0 2 Cl2 (SCOC) 
[J = OAt, t' = -0.34i, t" = 0.23i] and ft = 0.2i. The 
energies of the kink are 50-80 meV [l5[ and those of ac- 
tive phonons are 40-70 meV. Hence, with typical hopping 
amplitude t w QAeV one finds 0.12 < fl/t < 0.2. In 
further calculations we use spin-wave representation [55I ] 
and take into account only phonon-phonon vertex correc- 
tions. It was proved that the above approximations are 
reliable for the above parameters 5^, 5g|. 

The polaron picture for the finite doped case ex- 
plains the dichotomy between nodal and antinodal points 
[571 l58j | . The fact that the lifetime of the lowest peak in 
ARPES is large at the nodal point and small for the 
antinodal one, is a consequence of the contrasting na- 
ture of quasi-particle broadening by EPI at different en- 
ergies and momenta in the intermediate and weak cou- 
pling regimes. If the energy of the hole at momentum k, 
measured from the ground state of the Mott insulator in 
the nodal point is smaller than the phonon frequency ft, 
the decay of the quasiparticle by phonon emission is for- 
bidden by energy conservation law leading to the sharp 
peak. On the other hand, when the energy is larger than 
the phonon frequency (antinodal points, as an example) 
the real decay processes by phonon emission are allowed, 
causing significant line broadening even at moderate EPI 
coupling constants. For example, as it is seen from the 
calculated momentum dependence of the Lehman func- 
tion (Fig. Qli,b), the linewidth abruptly increases when 
the quasi-particle dispersion crosses the phonon energy 
at [k — tt/2]/tt rj 0.17. We emphasize that the interaction 
of holes with phonons is essential for the explanation of 
experimentally found "dichotomy". Although a broad- 
ening at general momentum k is present even in pure 
tt't" — J model, it is considerably smaller than the exper- 
imental one and, more important, does not show a sharp 



threshold-like increase above ~ 50 — 70meV. Of course 
one can not rule out inhomogeneous source of broaden- 
ing of antinodal quasiparticles (see, e.g. [Hil]) though the 
phonon mechanism is one of possible candidates to ex- 
plain the "dichotomy" . 

The low-energy dispersion (close to Fermi level sp), 
where quasiparticles are sharp, is separated from the 
high-energy part by a kink, where the velocity of the 
quasiparticle abruptly changes 15,24,[6(|. Experiments 
show that the high-energy velocity above the bosonic 
mode ft exhibits a strong doping dependence, while the 
doping dependence of the low-energy velocity is weaker 



24 1 . Fig. [2^, shows the experimentally extracted disper- 



sions from LSCO [24j]. These data were obtained using 



the momentum distribution curve (MDC) method, giv- 
ing the so-called universal low-energy velocity behavior 
where the velocity below the kink energy is found to be 
largely insensitive to doping. In the framework of the 
MDC analysis one traces the momentum dependence of 
intensity at a fixed energy. Then, the point of maximal 
intensity is considered to be an energy-momentum point 
representing the dispersion of the quasiparticle. 

DMC method is capable of obtaining the energy distri- 
bution curves (EDC) when one traces the energy depen- 
dence of the intensity keeping the momentum fixed. The 
reason is that the theoretical spectra here are obtained as 
analytic continuation to the real frequencies of imaginary 
time Green functions calculated by the DMC method at 
fixed momenta. Green functions calculated at different 
momenta possess slightly different distributions of statis- 
tical errorbars over the imaginary time. This difference 
is considerably amplified by the very high sensitivity of 
analytic continuation method to the distribution of sta- 
tistical errors. Hence, in contrast to EDC method the 
MDC method is not reliable in the framework of the the- 
oretical methodology used here. 

In order to compare directly experimental results with 
the calculated EDC dispersion, we also re-analyze the low 
energy data using the traditional EDC method. The ex- 
tracted MDC and EDC dispersions are plotted in Fig.[2ji 
and the corresponding low- and high-energy velocities are 
summarized in Fig. [3p. The MDC and EDC dispersion 
curves are very similar in the weak coupling limit, except 
for energies which are very close to phonon frequency [7| ■ 
However, MDC and EDC curves in the strongly interact- 
ing systems, as it was already observed before for bi-layer 
manganite [61|, are very different. The difference is es- 
pecially large at low energies while the high energy data 
are rather insensitive to the type of analysis. 

According to the analysis of the optical absorption (4|| , 
the EPI is enhanced with underdoping. Then, as is ex- 
pected in the Migdal-Eliashberg picture, the velocity at 
low energy should get smaller with underdoping while 
the high energy "bare" velocity above the phonon fre- 
quency should not change. The observed behavior is 
different. In the framework of the experimental MDC 
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FIG. l: (a) Evolution of the Lehman function along the nodal (tt/2, tt/2) — > (0, 0) direction for A = 0.5. 
Thick red lines indicate maxima of the Lehman function obtained from EDC analysis, (b) Intensity 
plot for the same parameters. Energy in panels (a-b) is counted from the vacuum state of the 
Hamiltonian (JTH2J) . (c) Momentum dependence of the peaks in the Lehman function for A = 0.126 
(circles), A = 0.32 (squares), and A = 0.50 (triangles). Energy in panel (c) is counted from the top of 
polaron band at k^v = (tt/2, 7r/2). The red arrow in (b) and (c) indicates the position of the phonon 
measured from the top of the polaron band. All energies are in the units of t — 0.4eV. 



analysis (Fig. [5b), similar to the theoretical EDC re- 
sults (Fig. [It), the high-energy velocity increases with 
decreasing doping (Fig. [5}d). While understanding the 
differences between EDC and MDC is a matter to be 
further explored, we compared EDC experimental result 
with EDC theoretical data which are only ones available 
from calculations done by DMC method. 

In contrast to the Migdal-Eliashberg picture, a polaron 
picture in the intermediate coupling regime can give a 
consistent explanation of these unique features. We cal- 
culated ARPES spectra for different values of A along 
the (0,0) - (n, 7r) nodal direction (Fig. QJ). Indeed, the 
energy of the ground state E(n/2, tt/2) of the Hamiltoni- 
ans (TJ312]) depends on the coupling constant A (Fig. [T^,,b) 
and, thus, the energies of the kinks at different As do not 
match to each other in the theoretical calculations. The 
energy onset in experiment is set in a different way. In the 
experimental data, presented in Ref. [24j , the binding en- 
ergy is counted from the Fermi energy ef and the origin 
of momentum is set to the Fermi momentum hp where 
the quasi-particle dispersion crosses the Fermi energy. In 



this manner the energy of the kink in experimental data is 
roughly the same for all dopings [24]]. Hence, to compare 
the results of our theoretical calculations with the exper- 
imental data for each A in Fig.[]Tc) we applied a constant 
energy shift to the theoretical results in order to match 
the energies of the quasi-particle for momenta just at the 
kink position. It is clear that the theoretical high-energy 
V ft f^ c (A) and low-energy C (A) velocities estimated 
just above and below the kink position, respectively, do 
not depend on the theoretical energy onset. However, 
both the high- and low-energy EDC velocities are depen- 
dent on coupling constant A, representing the intrinsic 
quasiparticle behavior. 

The calculated high-energy velocity (Fig. \Vp) in terms 
of the polaron picture gives a qualitative explanation of 
the anomalous doping dependence of the high-energy ve- 
locity observed in LSCO provided one assumes that the 
effective coupling constant A increases with underdop- 
ing. Indeed, one can see that, in contrast to the Migdal- 
Eliashberg picture, the theoretical high energy velocity 
Vfvigtfity (Fig-Hfc) just above the kink position increases 
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FIG. 2: (a). Energy- momentum dispersions for LSCO with different dopings, using both EDC and 
MDC methods. The MDC low (high) energy velocity V\ ow (V high ) is obtained by fitting MDC 
dispersion at binding energy 0~50meV (100~250meV) using a linear line. (b). Low and high- 
energy velocities as a function of doping obtained from MDC and EDC dispersions, (c). The EPI 
coupling constant obtained from EDC low-energy velocity and high energy velocity by empirical 
scaling relation © A = ^/(V^ h c (X) - V^ c '(A)) /V^ c (A) / V20 (filled circles). Empty squares and 

diamonds are coupling constants A 1 - 00 - 1 obtained from the analysis of the optical conductivity [45[ 
on weakly doped LSCO and YBa 2 Cu30 6+:! ; (YBCO), respectively. Only the ratio of A (oc) (x) at 
finite doping x and of that A (oc ^(x = 0) at zero doping is determined in Ref. [45j and we take 
A (oc) (a; = 0) = 0.85 to match results obtained from optical and ARPES measurements as close as 
possible. 



with A. Hence, according to experiment, the theoretical 
high energy velocity V^f^ (A) increases with underdop- 
ing. 

Furthermore, from our theoretical results of the spec- 
tral function of a single hole in the tt't" — J — ph 
model (TlJ[2j) , we can fit the dependence of the ratio 



LDC W - ^ C (A))/^l DC (A) on the effective con- 



(Vhigh 



pling constant A by a quadratic empirical scaling rela- 
tion: 



yEDC 
Low 



(A)) 



2Q\ Z 



(3) 



The doping dependence of the effective EPI constant A 
of LSCO estimated from the comparison of the above 
theoretical relation with experimental data is shown in 
Fig. [21;. Note, that the values of A obtained from ARPES 
are consistent with data derived from the doping depen- 
dence of the optical absorption 45] and other methods 



which take into account the strong interplay between EEI 
and EPI 0, [|| . We checked that the basic characteristics 
of the kink are mostly defined not by the range of EPI 
but by the resonance between the energy of the phonon 
and quasiparticle. Therefore, the data in Fig. [5£ are also 
valid for the Frohlich type of coupling. We note that 
our parameters are always in the large-polaron regime 
because EPI is less than the critical coupling A c = 0.6 
for given model [62| . The restriction of our analysis to 
a single phonon mode does not influence the relation Q 
because fine features, caused by interaction with m ultip le 
modes, do not change the gross shape of the kink [63| . 

Note here that the obtained effective As are reasonably 
of the order or less than unity even in the strongly under- 
doped regime where it is maximal. This is in sharp con- 
trast to the naive application of the Migdal-Eliashberg 
relation {Vhigh — Viow) /Viow = A which would give an 
unphysically large number on the order of ~ 10 at low 
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dopings. It is clear that the Migdal-Eliashberg approach 
underestimates the effects of EPI because the vertex cor- 
rections are neglected which, in turn, leads to an overes- 
timated value of the effective A. Indeed, even generalized 
Migdal-Eliashberg approach, including electron-electron 



correlation effects [6J, |65[, predicts that A ps 1.2 in the 



optimal doping, which is 3 times larger than our result 
obtained with full vertex corrections included. 

The values of A at low dopings, obtained in this paper 
from the kink angle, coincide with that obtained from the 
linewidth and distance of Franck-Condon peak from the 
chemical potential in undoped La 2 Cu04 compound [38j | . 



According to Fig. If in [15J, the kink angle is a rather 



universal function and, hence, the performed analysis 
applies also to La2-j / - x Nd 3/ Sr :c C04, Bi 2 Sr 2 CuOg and 
Bi 2 Sr 2 CaCu 2 8 . 

In the vicinity of the kink (Fig. QJi) high resolution 
Lehman spectral functions consist of the ground state 
peak and several phonon sidebands. The origin of the 
kink is the abrupt transfer of spectral weight from the 
low energy peak to the phonon sidebands when the mo- 
mentum is increased (Fig. QJi). The rate of transfer with 
momentum increase strongly depends on the EPI cou- 
pling constant and, thus, the velocity just above the kink 
energy is very sensitive to the EPI. 

Note, that the high energy velocity Vhigh is more sensi- 
tive to EPI change than the renormalized low energy ve- 
locity Vi ow (Fig.QJ). We note that the domain of the fast 
high energy velocity Vhigh is restricted to an energy range 
of a few phonon frequencies D, and to a small fraction of 
the Brillouin zone. Then, according to our calculations, 
at higher momenta/energies dispersion of the hump re- 
turns to the position of quasi-particle band noninteract- 
ing with phonons. The return of the quasi-particle dis- 
persion to the unrenormalized one is typical also for the 
weak coupling theory and for the Migdal-Eliashberg ap- 
proach neglecting the vertex corrections. However, in 
the intermediate coupling regime the energy/momentum 
domain of large velocity in exact DMC method is con- 
siderably wider than that obtained in a weak coupling 
approach. 

Several phonon sidebands survive above the momen- 
tum of kink (Fig.QJi) though they are not seen (Fig.Q})) if 
the spectrum is broadened more than it follows from the 
Hamiltonian tHH]) . Indeed there are numerous sources of 
decay, such as coming from hole-hole interaction, impuri- 
ties, surface roughness etc., which do not enter the Hamil- 
tonian (JT3C2J . The intensity map where numerical data 
are broadened by Gaussian with the width a — 0.1 does 
not show multiple sidebands but a single peak (Fig. [TJd) , 
which roughly corresponds to the energy of the highest 
intensity sideband. Comparing Figs. [T] and [2j the global 
features of the spectrum are well reproduced by our po- 
laron picture, which provides an alternative and com- 
plementary approximation replacing Migdal-Eliashberg 
approach in the underdoped regime. Strictly speaking, 



general features of the evolution of the Lehmann spectral 
function across the phonon energy are not restricted to 
the tt't" — J — ph model since similar phenomena should 
be observed for any polaronic model with intermediate 
strength of EPI [66]. What is novel here is that a similar 
scenario is realized in a more complicated situation with 
a combination of Mott physics and EPI. 

In conclusion, we propose a novel polaronic metal 
picture for underdoped cuprates and analyze the ARPES 
on underdoped cuprates in the framework of the pola- 
ronic scenario based on t-J model plus EPI. It is par- 
ticularly important that the electron-lattice polaron is 
formed in the background of the antiferromagnetic corre- 
lations and the essential influence of the electron-phonon 
interaction on the spectra is observed exclusively due 
to constructive interplay between electron-lattice and 
electron-electron interaction. Indeed, strong effects of 
coupling of holes to the lattice vibrations are observed 
only because the latter are highlighted by concomitant 
interaction of holes with short-range antiferromagnetic 
correlations. The latter ones persist up to moderate dop- 
ings even when the long-range antiferromagnetic order is 
destroyed. The doped holes form the small polaron in 
the undoped and heavily underdoped region. Analysis of 
ARPES in this scenario shows that, in accordance with 
the results previously obtained from the analysis of op- 
tical absorption, effective EPI decreases gradually and 
reaches the self-trapping crossover at optimal doping. 
The extended (large polaron) and localized (small po- 
laron) states coexist around this crossover region, which 
appear in a momentum dependent way. This picture 
explains many experimental puzzles, through the real- 
istic estimate of the coupling constant, such as the large 
momentum-dependent broadening of a single hole spec- 
tral function, the dichotomy of nodal and anti-nodal di- 
rection, and the unconventional doping dependence of 
ARPES of underdoped cuprates. 
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